Thermotropic Nematic Order Upon Nano-Capillary Filling 
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Optical birefringence and light absorption measurements reveal four regimes for the thermotropic 
behavior of a nematogen liquid (7CB) upon sequential filling of parallel-aligned capillaries of 12 nm 
diameter in a monolithic, mesoporous silica membrane. No molecular reorientation is observed for 
the first adsorbed monolayer. In the film-condensed state (up to 1 nm thickness) a weak, continuous 
paranematic-to-nematic (P-N) transition is found, which is shifted by 10 K below the discontinuous 
bulk transition at Tin =305 K. The capillary-condensed state exhibits a more pronounced, albeit 
still continuous P-N reordering, located 4 K below Tin. This shift vanishes abruptly on complete 
filling of the capillaries. It could originate in competing anchoring conditions at the free inner 
surfaces and at the pore walls or result from the 10 MPa tensile pressure release associated with the 
disappearance of concave menisci in the confined liquid upon complete filling. The study documents 
that the thermo-optical properties of nanoporous systems (or single nano-capillaries) can be tailored 
over a surprisingly wide range simply by variation of the the filling fraction with liquid crystals. 

PACS numbers: 61.30.Gd, 42.25.Lc, 64.70.Nd 



I. INTRODUCTION 

Liquid crystals spatially confined on the nanometer 
scale exhibit structural and thermodynamical proper- 
ties which differ markedly from their bulk counterparts. 
Both the collective orientational (isotropic-to-nematic) 
and translational (smectic) transitions have turned out 
to be significantly by finite size, quenched disorder and 
interfacial (solid-liquid) interactions introduced by the 
confining solid walls or the pore topology [l|-KLl|. 

For example, the first-order discontinuous isotropic-to- 
nematic transition of rod-like liquid crystals condensed 
in pores a few nanometers across goes to a critical point, 
above which, strictly speaking, it disappears. Then the 
transition is rendered continuous with a paranematic 
high-temperature phase, where even at highest temper- 
atures the pore walls (anchoring fields) induce a resid- 
ual collective orientational, paranematic order and no 
isotropic liquid is observable [J, 0, H, [l2| • 

This sensitivity of liquid crystalline phase transitions 
stimulated a plethora of experimental studies. Typically 
they were performed upon complete filling of the con- 
fining geometry [2, H, H-d) [l2| and thus revealed inte- 
gral quantities on the phase behavior. Phenomenologi- 
cal considerations [l| and computer simulations on liquid 
crystals in thin film and nanoporc confinement |13l4l8l | 
indicate, however, pronounced heterogeneities, in par- 
ticular interface-induced molecular layering and orienta- 
tional immobility in the pore wall proximity Despite 
recent experimental advancements in optical techniques 
directly probing orientational order parameter profiles in 
the proximity of planar, solid walls J19l42ll |. achieving the 
spatial resolution necessary to rigorously explore such in- 
homogeneities in nanometer-sized capillaries remains still 
experimentally extremely demanding. 



Here we present an experimental scheme which allows 
us to investigate the thermotropic orientational order of 
a classic, rod-like nematogen (7CB) as a function of grad- 
ual filling of parallel-aligned silica capillaries in a mono- 
lithic membrane. We find a gradual evolution of the ther- 
motropic orientational order behavior and can infer four 
distinct phase transition behaviors depending on the fill- 
ing fraction, and thus spatial arrangement of the liquid 
crystalline material in the capillaries. Moreover, we find 
hints of tensile pressure effects on the orientational be- 
havior of the nano-confined liquid crystal. 



II. EXPERIMENTAL 

In our experiments we employ a monolithic silica mem- 
brane of 320 ^m thickness. It is traversed by an ar- 
ray of parallel-aligned, non-interconnected capillaries. 
The membrane is prepared by thermal oxidation of a 
free-standing mesoporous silicon membrane [2 21 |23| at 
810°C for 12 h. The mean capillary diameter D and 
porosity P, as determined by recording a volumetric ni- 
trogen sorption isotherm at T =77 K, are 12.0 ±0.5 nm 
and 37 ±1 %, respectively 

In the bulk state 7CB exhibits an isotropic-to-nematic 
ordering transition at a temperature Tin=315 K. Given 
the low vapor pressure of 7CB at room temperature and 
the channel geometry with large, inner surfaces we can- 
not fill the matrix via the gas phase. Therefore, we im- 
bibe it with eight binary 7CB/acctone solutions [2J|. Af- 
ter evaporation of the high vapor pressure solvent, the 
matrix is filled with the remaining, low-vapor pressure 
liquid crystals up to a filling fraction /, that depends on 
the initial solute concentration. The filling fraction after 
evaporation is verified by weighting of the partially filled 
matrix and a comparison with the mass of the empty 



matrix of 20 mg. By a gradual increase of the solute 
concentration we prepare 10 filling fractions in the range 
from / =0.1 to 1 at room temperature. The evapora- 
tion of the solvent is controlled by time-dependent mass 
measurements of the sample at a temperature of 325 K 
and followed an exponential decay with time. After 1 
hour, the acetone had evaporated to 1/100 of its ini- 
tial concentration, and thus for the low-/ samples the 
further mass drop was below the detection limits of our 
balance, which was 5 ■ 10~^ g. Then we waited for an- 
other approx. 3 hours before we started the birefringence 
measurements in order to guarantee the absence of any 
significant amount of solvent in the sample. The esti- 
mated error margin in / due to the uncertainties in the 
mass determination amounted to 0.005. 



III. RESULTS AND DISCUSSION 

It is instructive to compare this filling procedure with 
the standard way for the preparation of partial liquid 
fillings of nanoporous substrates via the gas phase, i.e. 
vapor sorption isotherms. For nanoporous solids they 
exhibit distinct regimes typical of characteristic arrange- 
ments of the pore liquid: At small fractional filling /, 
below a certain threshold value /c, an adsorbed film is 
homogeneously formed at the pore walls, see sketch in 
Fig. [TJi, / < /c. For / > /c the pores fill via capil- 
lary condensation, that is the formation of liquid bridges 
terminated by concave menisci. Both regimes can be in- 
ferred by light adsorption experiments. Upon onset of 
capillary condensation or evaporation typically a dras- 
tic decrease in optical transmission, lg(///o) is observ- 
able |25l - [27| (Here Iq and / are the intensities of inci- 
dent and transmitted laser light, respectively). The ma- 
trix turns milky and the reduction in the transmitted 
light originates in scattering at the capillary bridges/void 
structures, which have sizes on the length scale of the 
wavelength of visible light [2^ [23]. By contrast, the ma- 
trix remains transparent both in the film condensed state 
(/ < /c) s-nd in the completely filled state (/ -^l), where 
such structures are absent. 

Laser light adsorption experiments on our solute- 
prepared samples exhibit an analogous behavior as a 
function of /, see Fig. [TJi. A sizeable decrease in light 
transmission is observed for / > /c = 0.25. The 
gradual disappearance of the bubbles and/or capillary 
bridges upon reaching a complete filling is indicated for 
/ > 0.95 — 0.97 and, eventually the full optical transmis- 
sion upon vanishing of all empty pore segments is recov- 
ered for / — ^ 1. Note that the transmission drops by four 
orders of magnitude in the capillary condensate regime, 
eventually even below the sensitivity of our setup. 

In contrast to nanoporous silica powders or disordered, 
sponge-like pore structures, our monolithic silica mem- 
brane is particularly suitable in order to evaluate the 
collective orientational order of liquid crystals by optical 
birefringence measurements in a transmission geometry 



[12l |. In the bulk nematic state the optical birefringence 
changes linearly with the orientational order parameter, 
Q = ^(3cos^ 6 — 1), where 6 is the angle between the 
long axis of a single molecule and a direction of preferred 
orientation of that axis, the director [30|. The brackets 
denote an averaging over all molecules under considera- 
tion. 

The simple relation between optical birefringence and 
Q the bulk system is also conserved for the liquid crys- 
tal/matrix composite material, provided the geometrical 
birefringence is subtracted [261 . Due to the changes in 
the geometrical arrangement of the liquid upon gradual 
filling, the latter one changes considerably in the entire 
region of partial fillings. This holds even in the case of 
the adsorption of a liquid with isotropic building blocks, 
like neopentane, which we experimentally and theoreti- 
cally demonstrated in two previous studies [2^ [3l|. It 
is, however, practically T-independent for a fixed /. For 
this reason the change of the optical birefringence, S{An) 
versus T, represents the quantity which most properly 
characterizes the orientational ordering inside the porous 
substrate. It was measured with a self-designed, high- 
resolution polarimetry setup presented in Refs. |12|.|32||. 
The high sensitivity of this device, which employs a dual 
lock-in detection scheme, allows one to obtain data even 
for low optical transmissions of the matrix. Therefore, it 
is possible to gain reliable data even upon onset of capil- 
lary condensation and upon reaching complete filling, as 
long as the transmission is larger than 10~^. 

In Fig. [2t S{An) is plotted versus temperature for a 
selected set of /'s in the film-condensed state. There 
exists a final birefringence and thus orientational order 
at high-T in all measurements. It is typical of a sizable 
collective ordering of the molecules induced by the in- 
teraction with the wall (interface anchoring). Silica sur- 
faces favor a planar anchoring in the case of the family 
of cyanobiphenyl liquid crystals (CBs), as documented 
in experiments j33l |34| and atomistic computer simula- 
tions [131. The magnitude and the positive sign of the 
changes of the optical birefringence for the lowest filling 
fraction (/ =0.12, corresponding to a partially filled, first 
adsorbed monolayer) indicates indeed such an alignment. 
It is even compatible with a preferred orientation of the 
long axis of the liquid crystals along the axial direction 
of the nano channels - as illustrated in Fig. [Hi. Note that 
there is a small, gradual increase S{An){T) upon cooling 
for this filling fraction, but no indication of a distinct 
change in collective orientational order. 

If one assumes a perfect cylindrical geometry and an in- 
termolecular distance of 0.4 nm 



35| 



the molecules added 
in the next filling steps correspond to ^-^1.5 and 2 mono- 
layers (/ =0.17 and 0.225). They exhibit sizeable in- 
creases in S{An){T) upon cooling, which are character- 
istic of a rearrangement between partial collective ori- 
entational order (a paranematic state) to full collective 
orientational order (nematic order). The symmetry of 
those states is identical; the transition is no "phase tran- 
sition" in a classical sense. It is, however, possible to dc- 
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FIG. 1: (Color online). 7CB in nanoporous silica, (a) Light transmission logio{I/Io) vs. filling fraction /, (b) Temperature 
of the paranematic-to-nematic transition TpN vs. filling fraction / as determined from the results of optical birefringence 
measurements shown in Fig. [21 (c) /-dependence of the effective optical birefringence (5(An)+ at T = TpN — 10 K as discussed 
in the text, (d) schematic side views of three characteristic capillary filling regimes: adsorbed monolayer(s), film regime (/ < /c), 
capillary condensate (/c < / < 1) and completely filled substrate (/ = 1). The molecular orientations in the different filling 
regimes are illustrated below (left panel d) and above (right panel d) TpN, that is in the nematic and paranematic state of the 
confined liquid crystal. In the panels (a-c) the /-range typical of capillary condensation is shaded. The solid lines in panel (a) 
and (b) are guide for the eyes, whereas in panel (c) linear fits are presented. 



fine TpN as a temperature of most pronounced changes 
of Q{T), in analogy to the "true" phase transition tem- 
perature. We determine this quantity by a calculation 
of the maximum in the T-derivative of S{An)(T) - see 
insets in Fig. [2j For the nanometric liquid crystal film 
it yields a TpN of 306 K. This downward-shift AT™ of 
lOK in respect to the bulk system compares surprisingly 
well with an expected ATjn of '^ 9 K, which one achieves 
by an extrapolation of thickness-dependent boundary ef- 
fects, explored in thin- film geometry on silica (with film 
thicknesses of several hundred nanometers) j36| . down 
to the nano metric films investigated here. The grad- 
ual increase in thcrmotropically induced molecular re- 
orientations upon film-thickness increase, observed here, 
corroborates a continuous increase in molecular orien- 
tational freedom with increasing distance from the pore 
walls documented in computer simulations for cylindrical 



nanopores [la, |l7[. 

The evolution of the thermotropic changes upon for- 
mation of capillary bridges (/c > / < 1) is documented 
in the panels (a) and (b) of Fig. [S] The transition from 
paranematic-to-nematic state remains continuous in this 
range of filling fractions whereas the temperature behav- 
ior of the optical birefringence and its relevant changes 
systematically scale with /. In order to characterize 
these changes quantitatively we introduce the quantity, 
(5(An)+ as the difference between the measured 6{An) at 
the lowest temperature investigated and its linearly ex- 
trapolated value from the paranematic phase far above 
Tpjy, see dashed line in Fig. [2b- By this subtraction we 
remove any contribution purely induced by the changes in 
the geometric birefringence upon capillary filling. There- 
fore, this quantity represents an effective order parame- 
ter Q = V~^ Jy Q{r)dV, characteristic of the integrated 



f: 1 - 0,120, 2 - 0.170, 3 - 0.225, 4 - 0.285, 
5-0.325, 6-0.385, 7-0.435, 8-0.925, 
9-0.950, 10-1.000 
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FIG. 2: (Color online). Temperature-dependent changes of 
the optical birefringence, 5{An), and its temperature deriva- 
tives (insets), —d{An)/dT, for selected filling fractions / as 
referenced in the header of the figure. Panels (a) and (b) 
correspond to the capillary-condensed regime and panel (c) 
represents the the film-condensed regime. Note that the wavi- 
ness of S{An) for f=0.925 results from the significant drop in 
transmitted light for that filling fraction - see Fig. 1. The 
dashed line in panel (b) examplifies a linear extrapolation of 
S{An) of the paranematic state, which is done for each filling 
fraction / in order to determine S{An)'^ . 



strength of the thermotropic, collective orientational (ne- 
matic) ordering inside the channels. 

In Fig. [T}: the /-dependence of 5{An)~^ is shown 
as determined in the confined nematic phase at T = 
Tin ~ lOK for each selected /. It increases linearly 
with / both in the film-condensed and in the capillary- 



condensed state. It is striking, however, that the slope in 
the experimentally accessible capillary-condensed regime 
(0.25<f<0.98) is by more than one order of magnitude 
larger than in the film-condensed state - see also in- 
set in Fig. [T};. This means that the film-condensed 
and capillary-condensed state is clearly distinguishable 
in terms of the evolution of thermotropic orientational 
order: In agreement with molecular simulation studies 
pTT ? ] , the freedom in orientational molecular ordering is 
considerably larger near the pore axis (i.e. in the core of 
the pore filling) than in the proximity of the channel wall, 
where the molecules are much more strongly influenced 
both by channel-wall roughness (quenched disorder) and 
the anchoring field of the channel walls |17l ? , |l3| • 

The linear increase of (S(An)+ with / in the film- 
condensed and capillary-condensed regime suggests that 
material added in one step does not affect the reorien- 
tational behavior of the material adsorbed in the pre- 
vious filling steps. We can not exclude, however, that 
the formation of capillary condensate in the pore cen- 
ter and the corresponding vanishing cylindrical liquid- 
vapor interfaces affects the reorientation behavior of the 
molecules in the film- condensed, boundary layer state. 

The temperature of the paranematic-to-nematic tran- 
sition exhibits only a weak increase in the entire region 
of capillary condensation (/c < / < 1) and remains at 
least ~ 4 K below the transition point TpN of the fully 
filled (/ = 1) porous substrate, see Fig. [TJd. 

Note, however, that there is a sizable upward jump of 
TpN between the slightly underfilled (/ < 1) and the com- 
pletely filled state (/ = 1) - see Fig. [TJ). Hence, Tpn of 
the continuous transition of the completely filled state is 
close to the temperature of the discontinuous bulk tran- 
sition Tin, in agreement with our previous experimental 
study and phenomenological treatment of this state with 
a Landau-de-Gennes model enriched by confinement and 
quenched disorder effects (KKLZ- model) jTI. la. [l3| . 

The abrupt T-shift in the transition temperature may 
result from the complete vanishing of any liquid/vapor 
interface and the resulting reduction of inhomogeneities 
in the system, which rather favor the disordered parane- 
matic state. In particular, for free liquid/vapor inter- 
faces, which are present for / < 1, the rod-like molecules 
are aligned parallel to the surface normal and the PN- 
transition in their vicinity is rather shifted upwards in T 
[37| . Albeit because of the high curvature of the menisci 
in the pores this anchoring condition competes with elas- 
tic energies, which effectively may lead to a reduction of 
TpN- Without having a detailed microscopic information 
of the orientational structure, a quantitative estimation 
of these influences appears difficult, in particular since at 
the triple line solid/vapor/liquid crystal of the menisci 
perimeters peculiar defect structures may occur. Molec- 
ular dynamics simulations on confined rod-like crystalline 
systems could be helpful to gain quantitative insights for 
this complex problem. 

There is another alternative explanation, which is re- 
lated to the hydrostatics of the confined liquid. In the 



capillary-condensed state the liquid experiences a ten- 
sile pressure, dictated by the concave curvature of the 
menisci terminating the liquid bridges. It is known that 
this negative hydrostatic pressure causes not only sub- 
tle deformations of the rigid, nanoporous matrix [38|-|40| , 
it also significantly affects density, and thus pressure- 
dependent first-order transitions, most prominently the 
liquid-solid transition [4ll443l |. Therefore, it is in princi- 
pal also expected to affect the isotropic-to-nematic tran- 
sition [4Jj Eal • According to the Laplace formula applied 
to the menisci in the capillaries, while assuming good 
wetting conditions (mean curvature radius = - pore ra- 
dius = -6 nm, surface tension of 7 CB=31.7 mN/m [4y|) 
the tensile pressure in the liquid bridges amounts to ~- 
10.5 MPa for /c < / < 1. This pressure is completely 
released upon reaching f — I (mean curvature radius ~ 
oo). An extrapolation of the pressure-dependence of Tin 
of 7CB, reported in the literature [ijl: towards negative 
pressures yields a corresponding downward shift of 3.8 K 
between partially and completely filled state, in good 
quantitative agreement with the observed T-shift. It is 
interesting to note that this T-shift is much larger than 
any tensile pressure effect on liquid crystals that has been 
reported to date. Presumably, this results from the ex- 
treme challenges ^^ associated with experiments aimed 
at reaching tensile pressures comparable to the ones es- 
tablished here in the partially filled state. Since the true 
curvature of the menisci at high /-values is not acces- 
sible in our experiments and since studies with higher 
f-resolution are experimentally demanding both due to 
the sizable light scattering in that regime and because 
of the error bars in the filling fraction (resulting from T- 
dependent volume changes in the liquid) , this conclusion 
remains somewhat speculative and more rigorous studies 
are planned for the future. 



IV. SUMMARY 

In summary our optical studies allowed us to experi- 
mentally document an inhomogeneous thermotropic ori- 



entational order behavior upon sequential filling of silica 
nano-capillaries with a rod-like liquid crystal which en- 
compasses four distinct regimes: the monolayer state, the 
film-condensed state, the capillary-condensed regime and 
the complete filling. 

Since oricntational ordering is key to functional materi- 
als with switching capability, we believe that our observa- 
tions are not only of high fundamental interest, but also 
of importance for the emerging field of nanophotonics 
[47| : Our study documents that the thermooptical prop- 
erties of nanoporous systems (or single nano-capillaries) 
can be tailored over a surprisingly wide range simply by 
variation of the the filling fraction with liquid crystals. 

For the future we envision analogous experimental 
studies on the molecular dynamics, which should show 
a similar complex behavior as the static properties inves- 
tigated here. We also hope that our findings stimulate 
temperature-dependent simulation studies and model 
calculations on partially filled nano-capillaries, aimed at 
an exploration of the remarkably complex phase transi- 
tion phenomenologies revealed here as a function of ma- 
terial distribution in the capillaries. 



ACKNOWLEDGMENT 



This work has been supported by the Polish Na- 
tional Science Centre under the Project "Molecular struc- 
ture and dynamics of liquid crystals based nanocompos- 
ites" (decision No DEC-2012/05/B/ST3/02782) and by 
the German Research Foundation (DFG) by Grant No. 
Hu850/3-l. 



[1] P. Sheng, Phys. Rev. Lett. 37, 1059 (1976). 

[2] G. Crawford and S. Zumber, eds., Liquid Crystals in 

Complex Geometries Formed by Polymer and Porous 

Networks (Taylor and Francis; New York, U.S.A., 1996). 
[3] T. Bellini, N. A. Clark, C. D. Muzny, L. Wu, X. Z., C. W. 

Garland, D. W. Schaefer, and B. J. Olivier, Phys. Rev. 

Lett. 69, 788 (1992). 
[4] G. S. lannacchione, G. P. Crawford, S. Zumer, J. W. 

Doane, and D. Finotello, Phys. Rev. Lett. 71, 2595 

(1993). 
[5] R. Guegan, D. Morineau, R. Lefort, A. Moreac, 

W. Beziel, M. Guendouz, J. M. Zanotti, and B. Frick, 

J. Chem. Phys. 126, 064902 (2007). 
[6] S. Qian, G. S. lannacchione, and D. Finotello, Phys. Rev. 

E 57, 4305 (1998). 



[7] Z. Kutnjak, S. Kralj, G. Lahajnar, and S. Zumer, Phys. 

Rev. E 68, 021705 (2003). 
[8] Z. Kutnjak, S. Kralj, G. Lahajnar, and S. Zumer, Phys. 

Rev. E 70, 051703 (2004). 
[9] C. Grigoriadis, H. Duran, M. Steinhart, M. Kappl, H. J. 

Butt, and G. Floudas, Acs Nano 5, 9208 (2011). 
[10] T. Araki, M. BuscagUa, T. Bellini, and H. Tanaka, Nat. 

Mat. 10, 303 (2011). 
[11] M. Ruths and B. Zappone, Langmuir 28, 8371 (2012). 
[12] A. V. Kityk, M. Wolff, K. Knorr, D. Morineau, R. Lefort, 

and P. Huber, Phys. Rev. Lett. 101, 187801 (2008). 
[13] T. Gruhn and M. Schoen, Phys. Rev. E 55, 2861 (1997). 
[14] T. Gruhn and M. Schoen, J. Chem. Phys. 108, 9124 

(1998). 
[15] K. Binder, J. Horbach, R. Vink, and A. De Virgihis, Soft 



Matter 4, 1555 (2008). 
[16] Q. Ji, R. Lefort, R. Busselez, and D. Morineau, J. Chem. [32 

Phys. 130, 234501 (2009). 
[17] Q. Ji, R. Lefort, and D. Morineau, Chem. Phys. Lett. [33 

478, 161 (2009). 
[18] A. Pizzirusso, R. Berardi, L. Muccioli, M. Ricci, and [34 

C. Zannoni, Chem. Science 3, 573 (2012). 
[19] V. Barna, A. De Luca, and C. Rosenblatt, Nanotechnol- [35 

ogy 19, 325709 (2008). 
[20] A. De Luca, V. Barna, T. J. Atherton, G. Carbone, M. E. [36 

Sousa, and C. Rosenblatt, Nat. Phys. 4, 869 (2008). 
[21] J.-H. Lee, T. J. Atherton, V. Barna, A. De Luca, [37 

E. Bruno, R. G. Petschek, and C. Rosenblatt, Phys. Rev. 

Lett. 102, 167801 (2009). [38 

[22] S. Gruener and P. Huber, Phys. Rev. Lett. 100, 064502 

(2008). [39 

[23] P. Kumar, G. Franzese, and H. E. Stanley, Phys. Rev. 

Lett. 100, 105701 (2008). [40 

[24] S. Gruener and P. Huber, J. Phys. : Cond. Matt. 23, 

184109 (2011). [41 

[25] J. H. Page, J. Liu, B. Abeles, H. W. Deckman, and D. A. 

Weitz, Phys. Rev. Lett. 71, 1216 (1993). [42 

[26] A. V. Kityk, K. Knorr, and P. Huber, Phys. Rev. B 80, 

035421 (2009). [43 

[27] V. P. Soprunyuk, D. Wallacher, P. Huber, K. Knorr, and 

A. V. Kityk, Phys. Rev. B 67, 144105 (2003). [44 

[28] S. Naumov, A. Khokhlov, R. ValiuUin, J. Kaerger, and 

P. A. Monson, Phys. Rev. E 78, 060601(R) (2008). [45 

[29] S. Naumov, R. ValiuUin, J. Kaerger, and P. A. Monson, 

Phys. Rev. E 80, 031607 (2009). [46 

[30] S. Kumar, ed., Liquid Crystals (Cambridge University 

Press, Cambridge (UK), 2001). [47 

[31] M. WolflF, K. Knorr, P. Huber, and A. V. Kityk, Physical 



Review B 82, 235404 (2010). 

M. Skarabot, M. Cepic, B. Zeks, R. Blinc, G. Heppke, 

A. V. Kityk, and I. Musevic, Phys. Rev. E 58, 575 (1998). 

F. Vandenbrouck, S. Bardon, M. P. VaUgnat, and A. M. 
Cazabat, Phys. Rev. Lett. 81, 610 (1998). 

L Drevensek Olenik, K. Kocevar, L Musevic, and T. Ras- 
ing, Eur. Phys. J. E 11, 169 (2003). 

A. V. Komolkin, A. Laaksonen, and A. Maliniak, Journal 
of Chemical Physics 101, 4103 (1994). 
H. Yokoyama, J. Chem. Soc. - Faraday Trans. li 84, 1023 
(1988). 

H. Kasten and G. Strobl, Journal of Chemical Physics 
103, 6768 (1995). 

G. Guenther, J. Prass, O. Paris, and M. Schoen, Phys. 
Rev. Lett. 101, 086104 (2008). 

J. Prass, D. Mueter, P. Fratzl, and O. Paris, Appl. Phys. 

Lett. 95, 083121 (2009). 

G. Y. Cor and A. V. Neimark, Langmuir 26, 13021 

(2010). 

K. Morishige and H. Yasunaga, J. Phys. Chem. B 110, 

3864 (2006). 

C. Schaefer, T. Hofmann, D. Wallacher, P. Huber, and 

K. Knorr, Phys. Rev. Lett. 100, 175701 (2008). 

S. T. Moerz, K. Knorr, and P. Huber, Phys. Rev. B 85, 

075403 (2012). 

C. Reinecke and D. Demus, Cryst. Res. and Techn. 28, 

273 (1993). 

V. Manjuladevi, R. Pratibha, and N. V. Madhusudana, 

Phys. Rev. Lett. 88, 055701 (2002). 

U. Delabre, C. Richard, and A. M. Cazabat, J. Phys.: 

Cond. Matt. 21, 464129 (2009). 

I. Abdulhalim, J. Nanophotonics 6, 061001 (2012). 



